Exact relations between particle fluctuations and entanglement in Fermi gases 
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We derive exact relations between the Renyi entanglement entropies and the particle number 
fluctuations of (connected and disjoint) spatial regions in systems of A'^ noninteracting fermions in 
arbitrary dimension. We prove that the asymptotic large- A'^ behavior of the entanglement entropies 
is proportional to the variance of the particle number. We also consider ID Fermi gases with a 
localized impurity, where all particle cumulants contribute to the asymptotic large- A^ behavior of 
the entanglement entropies. The particle cumulant expansion turns out to be convergent for all 
integer-order Renyi entropies (except for the von Neumann entropy) and the flrst few cumulants 
provide already a good approximation. Since the particle cumulants are accessible to experiments, 
these relations may provide a measure of entanglement in these systems. 
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The nature of the quantum correlations of many- 
body systems, and in particular the entanglement phe- 
nomenon, are fundamental physical issues. They have at- 
tracted much theoretical interest in the last few decades, 
due to the impressive progress in the experimental activ- 
ity in atomic physics, quantum optics and nanoscience, 
which has provided a great opportunity to investigate 
the interplay between quantum and statistical behaviors 
in particle systems. The great ability in the manipula- 
tion of cold atoms in optical lattice (see, e.g., Ref. [l|) 
has allowed the realization of physical systems which are 
accurately described by theoretical models such as Hub- 
bard and Bose-Hubbard models in different dimensions, 
achieving through experimental checks of the fundamen- 
tal theoretical paradigma of condensed matter physics. 

The quantum correlations arising in the ground state 
of quantum many-body systems can be characterized by 
the expectation values of the products of local operators, 
such as the particle density and one-particle operators, 
or by their integral over a space region A, such as the 
particle-number correlators within A, 
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where n(x) is the particle-density operator and 
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counts the number of particles in A. Quantum corre- 
lations arc also characterized by the fundamental phe- 
nomenon of entanglement, which gives rise to nontrivial 
connections between different parts of extended quantum 
systems [1] . A widely accepted measure of entanglement 
is given by the Renyi entropies of the reduced density 
matrix pA of a subsystem A: 
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whose limit a —J- 1 provides the von Neumann (vN) en- 
tropy. Local correlations and bipartite entanglement en- 



tropies provide important and complementary informa- 
tion of the quantum features of many-body systems, of 
their ground states and of their unitary evolutions, be- 
cause they probe different features of the quantum dy- 
namics. However, the entanglement entropy is a highly 
nonlocal quantity which is difficult to measure. Designing 
an experimental protocol for its measurement represents 
a major challenge. 

A recent interesting proposal considers the particle 
fluctuations as effective probes of many-body entangle- 
ment This is based on the result that, for non- 
interacting fermions, one can write down a formal ex- 
pansion of the entanglement entropies of bipartitions in 

(2k) 

terms of the even cumulants of the particle-number 

distribution, which can be defined through a generator 
function as 

y(")==(-z5;,)"ln(e*^^-^)U=o. (4) 
Indeed; the Renyi entropies can be written as [3, 0] 
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where C is the generalized Riemann zeta functioir. In 
particular, for the lowest integer a we have 

- ^A'' + ^T^.vf + - (7) 
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The possibility of turning these expansions into an effec- 
tive measure of entanglement depends on its convergence 
properties, which appear problematic due to the behavior 
of the coefficients s^"'' with increasing k. 

In this paper we investigate the relations between en- 
tanglement entropies and particle fluctuations, and in 
particular the convergence properties of the formal ex- 
pansions ([S][5|). We will show that, in noninteracting 
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fermion gases with N particles in a finite volume of any 
dimension d, the expansion ([S]) gets effectively truncated 

(n) 

in the large-A^ limit, because the high cumulants 
with n > 2 arc all suppressed relatively to the particle 
variance 
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The leading TV''' ^^^'' inN asymptotic behavior of S*^"-* 
in Eqs. JS]), d?]) and dS]) arises from ' only, because 
the leading order of each cumulant V^'^'' with k > 2 van- 
ishes for any subsystem A (including disjoint ones) in any 
dimension. This implies the general asymptotic relation 



0(1). 
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We will also consider the effect of a localized impurity 
in ID fermion gases. In this case all particle-number 
cumulants contribute to the leading logarithmic behavior 
of the Renyi entropies. The expansion ([5]) turns out to 
be convergent for the integer Renyi entropies, so that a 
few particle cumulants provide a good estimate, with the 
exception of the vN entropy which appears problematic 
in this respect. Relations like Eq. (jlOp may provide an 
experimental measure of entanglement. 

This letter is organized as follows. First we consider 
Fermi gases in any dimension, and prove Eq. (|10p for 
any subsystem A, by rigorously computing the asymp- 
totic behaviors of all cumulants of the particle distribu- 
tion. Then we focus on one-dimensional Fermi gases, for 
which more general results, including subleading terms, 
can be obtained for the asymptotic behaviors of the parti- 
cle fluctuations and entanglement entropies of connected 
and disjoint subsystems. Finally, we extend our analy- 
sis to one-dimensional Fermi gases in the presence of a 
defect. 

The leading behavior in arbitrary dimension. - 

We consider a system of N non-interacting spinless 
fcrmions with discrete one-particle energy spectrum, 
which may arise from a finite volume or an external po- 
tential. The many-body ground state is obtained by fill- 
ing the lowest N one-particle energy levels (/)„(a;). The 
cumulants V^'^ of the particle-number distribution and 
the entanglement entropies of a subs;^stem A can be writ 
ten in terms of the overlap matrix 



syster 

m 



d'*z</.;(z)0™(z), n,m = l,...,iV. (11) 



The eigenvalues of A are real and limited, G (0, 1). 
The matrix 



I- A 
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is the overlap matrix of the complement of the region A. 

The particle fluctuations within a region A can be char- 
acterized by the cumulants 1/^'"' of the particle distribu- 
tion. The cumulant generator function for lattice free 



fcrmions has been derived in Ref. |6| in terms of the two- 
point correlation function Ca restricted to the subsystem 
A. From this, taking the continuum limit and using the 
fundamental property 



TtCX = TrA" 



for any n G N, 
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the cumulants generator function of the particle distri- 
bution can be expressed in terms of the overlap matrix 
A by replacing Ca with A, obtaining 



-(-*aA)"G(A,A)u^o, 

G(A,A) =Trln[l+(e'^-l). 
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The even cumulants which enter Eq. ([5]), can be 

cast in the form 



where 
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E = A(l - A) = AA, 



In particular, 

Vj^^ = TrE, V[^'^ = Tr [E - eE^]. (19) 
The entanglement entropies are obtained as 0, H] 

5^") = -J-_Trln[A" + (1 - A)"]. (20) 
In particular, the a = 2 entropy can be written as 
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The eigenvalues of E satisfies G (0, 1/4), thus the se- 
ries (|2ip is convergent for any iV, providing a system- 
atic approximation scheme in terms of vj^'^'^ by inverting 
Eq. 1111). 

In systems of noninteracting fcrmions and for arbitrary 
dimension d, the entanglement entropy of connected bi- 
partitions grows asymptotically as N^'^~^^/'^ In 0. The 
logarithm of the asymptotic behavior is related to the log- 
arithmic area-law violation in lattice free fermions |10l - 
[l6| . In homogeneous systems with periodic and open 
(hard-wall) boundary conditions (PBC and OBC respec- 
tively) the prefactor can be analytically computed [9| us- 
ing the Widom conjecture [l^l- This method applies, and 
is even more suited, to compute the large- behavior of 
particle fluctuations with both PBC and OBC in any di- 
mension. In fact, unlike the entanglement entropies, we 
deal with smooth functions, for which the Widom con- 
jecture has been proved [l^, [l^. We apply this theorem 
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to the overlap matrix of a subsystem A (with smooth 
bomidaries dA) of a finite system of size L'' with PBC, 
which is [§] 

A„„ = L-"^ [ df^x e»2'^(kn.-k„).x/L (22) 

J A 

with k G Z'' within the Fermi surface dT. It ahows us 
to derive the large- iV behavior of TrF(A), where F{z) is 
any function analytic in {z : |z| < 1 + e} with F(0) = 
F{1) = 0, obtaining 
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where we set L = 1, rix and are the normal vectors 
on dA and on df which is the Fermi surface dT rescaled 
to enclose a unit volume, and 
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Note that the function F enters only the integral I{F). 
The result ^ applies also to OBC. 

The asymptotic large- iV behavior of Tr E*"' can be com- 
puted using Eqs. (|23p and (|24|) . Indeed, it corresponds 
to the function 



Ffe(z) = z'=(l-z)* 



thus 
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Plugging the last result into Eq. (|23p , we finally have 



TrE^ =7Vi-hn7V- ^^'1' 
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The large- iV leading behaviors of V^^ are obtained by 
inserting these asymptotic results into their expressions 
in terms of TrE*^. cf. Eq. ((T6)) . For any spatial region A 

(2) 

in any dimension d, the variance V\ is 
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while, very remarkably, this leading term cancels for 
higher cumulants. For odd cumulants y(2fc+i) fj^g leading 
term vanishes, because they arc odd under A — >■ I — A. 

In the sum ([5]), the leading behavior of the Renyi en- 
tropies gets a finite contribution only from the variance 
and the resulting entropies agree the direct computation 
in Ref. @. Taking the ratio, the asymptotic large- rela- 
tion ([TU)) follows. The above calculations do not allow us 
to determine the behavior of the suppressed corrections 



in Eq. dTO]). Finite- iV calculations up to TV = O(IO^) in- 
dicate that they are 0(1/ In A^) . This is also supported by 
the analytic calculations in ID systems reported below. 

Notice that Eq. ([TU]) relating particle fluctuations and 
entanglement entropies can also be obtained for lattice 
free fcrmions in the thermodynamic limit, exploiting the 
correspondence [1, 0| between the overlap matrix A and 
the lattice two-point function where i, j are the lattice 
sites within the region A. Indeed, analogous formulas for 
the particle-number cumulants hold by replacing A with 
C and the lattice version of some of the above results 
has been already derived (as e.g. in Refs. M [ill , [isj). 
Moreover, developing the results of Refs. [1^, [201, anal- 
ogous results can be also inferred for free fermion gases 
in external potential, such as an harmonic one which is 
usually present in experiments of cold atoms [l|. 

One-dimensional Fermi gas.- For ID systems one 
may also consider an alternative computation based 
on the Fisher-Hartwig conjecture [2l| and generaliza- 
tions [2^ , similarly to what has been done for the en- 
tanglement entropies H, [H, [2^. This exact approach 
allows us to calculate not only the leading term in V^™-* 
but also the leading and subleading corrections to the 
scaling. 

Let us consider a ID system of size L with PBC or 
OBC, and the interval A = [0,x] as subsystem. Setting 
L = 1, the corresponding overlap matrix is 01 



sin[7r(n — m)x] 
7r(n — m) 



(PBC) (29) 



and 



sin[7r(n — m)x] sin[7r(ri + 771)2;] 



7r(n — m) 



7r(n -|- m) 



(30) 
(OBC). 



The asymptotic large-A^ behavior of the entanglement 
entropies has been already reported in Refs. 0, [E] 



s^'^Hx) = 

for PBC, and 

S^°''>{x) = 



a 



■ln(iVsin7rx) H-fea + 0(1) (31) 
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■ ln(2iVsin7rx) 



+ o(l) (32) 



for OBC; where also the constants 6^ and the o(l) cor- 
rections are known. Fisher-Hartwig calculations can be 
generalized to the particle cumulants (details will be re- 
ported elsewhere), obtaining for PBC 



TrE'' 



with 



1 [ik-l)l? 
7r2 (2/0-1)! 



ln(27Vsin nx) 



(33) 



5fe(tanW) r(l/2 + zx) 

+ il dx -J- — -In—— — -Fol, 

cosh (ttx) r(l/2-za;) 

gk{z)^kz[{l-z^)/Af-^/2. (34) 
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The leading terms agree with Eq. (|27)) , but Eq. (|33)) also 
provides the 0(1) corrections. Then, using Eq. (jl6p . we 
obtain 

= TrE = ^ln(A^sin7r2;) + V2 + 0(1), (35) 



V 



where 



[x) = V2k + 0(1) for k > 1, 



V2 



(l+7£ + ln2)/^2 
-0.0185104, U6 - 



0.00808937, 
etc. The odd cumulants are suppressed, 

y(2fc+l)(2;) =.0(1). 

For OBC we find 

V^'^l{N)^vj'Sli2N)/2 + o{l). 



(36) 



(37) 
(38) 



(39) 



(40) 



Following Ref. [2J], the o(l) subleading corrections can 
be systematically obtained. 

Asymptotic behaviors for disjoint subsystems. 

We now show that the asymptotic relation (fTO|) holds 
also for disjoint subsystems, for which some exact results 
can be obtained. Let us consider a Fermi gas of size 
L = 1 with PBC. The overlap matrix of a generic disjoint 
subsystem 

[a;i,a;2] U [2:3, 2:4], < xi < .T2 < < a;4 < 1, (41) 
reads 

A„™ = e'^("-'")(^^+^-^)p„™(x4 - xa) (42) 
+ e"("-™)(^^+^i)p„„,(x2-xi). 



where the matrix P„„i is defined in Eq. ([29]). For the 
particular subsystem 

B= [0,a;]U[l/2,l/2 + a;], < x < 1/2, (43) 

the overlap matrix takes the simple form 

A™„ = e-("^")^[l + (-l)"-"]P,„„(x). (44) 

Since we are interested in the traces of powers of A„m , we 
may neglect the global phase. Then, exploiting the fact 
that all terms with different parity (odd n — m) vanish, 
we have for any integer k 



Tr 



,(iV) = 2TrA^(iV/2), A = [0,2a;], 



(45) 



where Aa{N/2) is the overlap matrix of the connected in- 
terval A = [0, 2x] with N/2 particles reported in Eq. (f29| . 
This exact relation implies 

S^s\N) = 2S^f{N/2), (46) 
V^"'\n) = 2Vj-"\N/2). (47) 



Using the asymptotic behaviors ([3T|). ([35]) and (|36l) for 
the r.h.s. of the above equations, we again recover 



the asymptotic relation (|10p . Moreover, we obtain the 
asymptotic behaviors of the entanglement entropies 



1 + a 



■ (lniV + lnsin27rx) +6a + 0(1), (48) 



which agree with the free-fermion conformal field theory 
prediction [2^ 



-1 r- 
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[i:i,a:2]U[2:3,a;4] ~ g 

sin(7r.T2i) sin(7ra;43) sin(7ra;4i) sin(7ra:32) 



-In: 



+6a + 0(1), 



sin(7rx3i) sin(7ra;42) 



(49) 



where Xij = Xi — xj and < xi < 0:2 < a;3 < 2:4 < 1. 
While we explicitly reported the proof of Eq. ^TU\\ for the 
particular choice of the intervals in Eq. the proof 

can be extended to subsystems with an arbitrary union 
of equal and equidistant intervals. We expect that the 
validity of Eq. ([TU]) extends to the most general case of 
disjoint subsystems. 

One-dimensional Fermi gas with an impurity.- 
We now investigate how the above relations between en- 
tanglement entropy and particle fluctuations may change 
in the presence of localized interactions, such as those 
arising from local impurities. We consider ID fermion 
gases with an impurity localized at the point separating 
the system in two equal parts of size L ~ 1, with hard- 
wall boundary conditions at their ends. The allowed scale 
invariant conditions at the vertex, describing the univer- 
sal features arising from the presence of a defect, are fully 
encoded in the scattering matrix [2^ 



m = 



l + e^ 



-l + e^ 2e 
2e 1 - e2 



T = 



2e 



l + e^ 



(50) 



where T is the transmission coefficient. T = I corre- 
sponds to full transmission, i.e., no impurity, thus to 
the bipartition into two equal parts of a ID system with 
OBC. The corresponding entropies and cumulants can be 
obtained by setting a; = 1/2 in Eqs. p2j|36p . 

For \T\ < 1 the ground-state entanglement entropies 
and particle fluctuations of one of the two edges can be 
again derived from its overlap matrix A. For even N we 
have 0, [13 



Km = , ^'^ 0„m(l/2) for n 7^ m, 
1 + e'^ 

1 £2 

^ - — — 5- for odd 71, A„„ = — — - 

1 + 1 + 



(51) 



for even n. 



where 0„m is defined in Eq. pop . The symmetry ^ 
1 — Ofc of the spectrum of A implies that any odd observ- 
able with respect to A ^ 1 — A, such as the odd particle 
cumulants, vanishes. 

Using the result of Ref. [13 , 



TvWiT) = r'"'TrE''(T = 1) 



(52) 
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TABLE I: In order to check the convergence of the expan- 
sions (O for the coefficients Cg(a) of the large-Af logarith- 
mic behavior S^'^'> « C^t^) in TV in the presence of a defect 

with T = 1/2, we report the sum •Sfc^^C'vf^'') f^o'" ^^'^ 

a = 1, 2, 3 Renyi entropies, whose leading-log coefficients are 
0.0570281, 0.0264623, 0.0203204 respectively. An analogous 
pattern of convergence is shown by the data at fixed A''. 



K 


vN (a = 1) 


a = 2 


a = 3 


1 


0.0416667 


0.0312500 


0.0277778 


2 


0.0622283 


0.0264309 


0.0201623 


3 


0.0622283 


0.0264309 


0.0201623 


4 


0.0288468 


0.0264614 


0.0203150 


5 


0.0763883 


0.0264626 


0.0203256 


6 


0.703467 


0.0264623 


0.0203209 


7 


-3.61878 


0.0264623 


0.0203202 


8 


-47.3035 


0.0264623 


0.0203203 


9 


949.44 


0.0264623 


0.0203204 


10 


6860.34 


0.0264623 


0.0203204 



the particle fluctuations and entanglement entropies can 
be computed using Eqs. (jl4j|21[) . In particular, for the 
particle variance we obtain 

y(2)(y) = TrE(r) = T2y(2)(T = 1). (53) 

This shows that, like the homogeneous case T = 1, the 
particle variance grows as In . but with a smaller coeffi- 
cient r^/ (27r^). Analogous results for the large- iV behav- 
ior of higher cumulants follow and, unlike the homoge- 
neous case, also the higher-order even particle cumulants 
grow logarithmically when |r| < 1 

Cy(2.,(T)lnAf. (54) 

Therefore, the relation ^ between entanglement en- 
tropies and cumulants does not get truncated at large N 
and Eq. PH]) is not valid anymore. However, the large 
order behavior in Eqs. and ([5]) show that these in- 
finite sums are convergent for any integer a > 1, while 
for non-integer a they are not convergent. The truncated 
sum of the leading large- A^ logarithms of the particle cu- 
mulants rapidly approaches the value of the coefficient of 
the integer Renyi entropy, see Table [H with the only ex- 
ception of the vN entropy. In particular, using Eq. (1211) . 
we obtain 

5(2) = C5(2,(T)lniV + 0(l), (55) 

The comparison with the entanglement entropy and 
particle variance of two hard-core Bose-Hubbard (XX) 
chains of size L separated by a bond defect [H, [29| . 
5(2) ~ Cs(2,(r)lnL and 1/(2) ^ Cyi2){T)\nL respec- 
tively, shows the universality of the leading logarithmic 
behaviors. 



Conclusions. - 

We showed that, in a noninteracting Fermi gas with a 
large number N of particles and in any dimension, the 
Renyi entanglement entropies of a spatial (connected or 
disjoint) region A are proportional to the variance of the 
particle number in the same region, cf. Eq. (jlO[) . with 
a coefficient that does not depend either on A or on the 
space dimension of the system. This remarkable result 
is proved by first computing the asymptotic behaviors 
of all particle cumulants V'^™\ which are such that the 
leading order of each cumulant V'^^^ with k > 2 van- 
ishes for any subsystem A (including disjoint ones) in 
any dimension. Then, the leading iV^''"^-'/'' IniV asymp- 
totic behavior of 5^"' in the formal expansion S'^"^ = 

J2kLi ^i"''^!^'^'' arises from V^'^ only. Therefore, en- 
tanglement entropies are directly related to the integral 
of the density-density correlation function which is ac- 
cessible in cold atoms experiments through Bragg spec- 
troscopy (30l |. thus providing an experimental measure- 
ment of the zero-temperature entanglement. Further- 
more the recent developed technique of single atom-site 
imaging technique [sH should allow to measure all the 
lowest particlc-numbcr cumulants. 

The situation is more involved for interacting systems. 
Already, in systems with localized interactions arising 
from impurities, all the cumulants T/t^'') contribute to 
the asymptotic large- A^ behavior of the entanglement en- 
tropies in the expansion S'^"'' = X]fc°=i ■^i"''^!^'^'' ■ How- 
ever, the resulting sum turns out to be rapidly conver- 
gent for all integer Renyi entropies (except for the vN 
one), and a few cumulants are enough to get a precise 
estimate of the entanglement. 

The conservation of a global charge, and in particular 
the particle number, is crucial for the connections be- 
tween bipartite entanglement and particle fluctuations. 
For interacting systems not conserving the particle num- 
ber, the entanglement should be related to the more fun- 
damental energy transport. Recent proposals d, HI, H^l 
to measure the entanglement entropy in general ID sys- 
tems have considered protocols based on appropriate 
quantum quenches. 

Finally, let us mention that the relations between parti- 
cle fluctuations and entanglement entropies apply also to 
ID Bose gases with short-ranged repulsive interactions. 
Indeed, in the limit of strong interaction (i.e. a gas of 
hard-core bosons), this can be mapped to free fermions, 
so that their particle-number cumulants and entangle- 
ment entropies coincide for connected regions. The hard- 
core Bose gas also describes the dilute limit of the flnite- 
strength models [s^], which implies that it also provides 
their infinite size limit at fixed N [2^ (with 0{L~^) cor- 
rections). Therefore, in this regime, the same asymptotic 
large- A^ behavior is expected, and in particular Eq. ([TU]) 
remains valid. 
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